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Abstract 

We develop a quasiclassical theory of the synchrotron-hke gluon radiation. Our calcu- 
lations show that the parton energy loss due to the synchrotron gluon emission may be 
important in the jet quenching phenomenon if the plasma instabilities generate a suffi- 
ciently strong chromomagnetic field. Our gluon spectrum disagrees with that obtained 
by Shuryak and Zahed within the Schwinger's proper time method. 

1. The hydrodynamic description [1] of the RHIC hadron spectra at low pt gives strong 
evidence for the thermalization of the produced hot quark-gluon plasma (QGP) at the 
time scale tq ~ 1 fm. The early thermalization provides a serious challenge for the per- 
turbative approach based on the Boltzmann equation with the binary and the 2^3 
processes [2]. The fast thermalization can probably be explained if the plasma is strongly 
coupled [3]. A promising possibility of speeding the equilibration in a weakly coupled 
plasma scenario is connected with plasma instabilities [4, 5, 6, 7, 8]. One of the poten- 
tially important instabilities is the Weibel [9] instability which occurs for the anisotropic 
plasma distribution. In the QGP produced in y4y4-collisions the initial parton distribution 
has Pz <^ Pt [2]. In this case the Weibel instability can generate considerable (predom- 
inantly transverse) chromomagnetic field which bends the parton trajectories and leads 
eventually to isotropization of the parton distribution [5, 7, 8]. Also, the generated chro- 
momagnetic field may be important for explanation of the anomalously small plasma 
viscosity [10] and the longitudinal broadening of the jet cone (the ridge effect) [11]. 

In the scenario with generation of the chromomagnetic field an interesting ques- 
tion arises on the effect of the chromomagnetic field on the jet quenching due to the 
synchrotron-like gluon emission from fast partons. Evidently the available analyses of the 
jet quenching, based on the idea of the induced gluon emission due to multiple scattering 
[12, 13, 14, 15, 16], become inappropriate if the contribution of this mechanism to the 
parton energy loss is comparable to that from the ordinary multiple scattering. It is clear 
that in this case both these mechanism must be treated on an even footing. As a first step 
in understanding the role of the chromomagnetic field in the parton energy loss it would 
be interesting to calculate the purely synchrotron induced gluon radiation neglecting the 
interference between the two mechanisms. It requires a formalism for the synchrotron 
radiation in QCD. 

In the present paper we develop a semiclassical approach to the synchrotron gluon 
radiation and give a qualitative estimate of the parton energy loss in the QGP produced 
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in 7474-collisions. For photon emission our result agrees with prediction of the quasiclassical 
operator approach [17]. For gluon emission our spectrum disagrees with that obtained 
by Shuryak and Zahed [18] within the Schwinger's proper time method. We give simple 
physical arguments that the spectrum obtained in [18] is incorrect. 

2. We consider the synchrotron gluon radiation from a fast parton in the quasiclassical 
regime when for each parton (initial or final) the wavelength is much smaller than its 
Larmor radius, Rl- One can show that in this regime, similarly to the photon radiation in 
QED, the coherence length of the gluon emission, L^., is small compared to the minimal Rl. 
It allows one to perform the calculation of the radiation rate per unit length by considering 
the case of a slab of chromomagnetic field of thickness L which is large compared to L^, 
but small compared to the minimal Rl. In this case the transverse momenta of the final 
partons are small compared to their longitudinal momenta (we choose the 2;-axis along 
the initial parton momentum) . We consider the case of a slab perpendicular to the z-axis 
with transverse chromomagnetic field, H^. It is enough to consider the chromomagnetic 
field with the only nonzero color components in the Cartan subalgebra, i.e., for a = 3 
and a = 8 for the SU (3) color group. The interaction of the gluons with the background 
chromomagnetic field is diagonalized by introducing the fields having definite color isospin, 
Qa, and color hypercharge, Qb, (we will describe the color charge by the two-dimensional 
vector Q = {Qa, Qb))- In terms of the usual gluon vector potential, G, the diagonal color 
gluon states read (the Lorentz indices are omitted) X = (Gi +iG2)/V^ {Q = (—1,0)), 
y = (6-4 + iG,)/V2 {Q = (-1/2, -x/3/2)), Z = {Ge + iG,)l^ {Q = (1/2, -%/3/2)). 
The neutral gluons A — G^ and B — G3 with Q = (0, 0), to leading order in the couphng 
constant, do not interact with the background field, and the emission of these gluons are 
similar to the photon radiation in QED. 

The S'-matrix element of the q gq' synchrotron transition can be written as (we 
omit the color factors and indices) 

{gq'\S\q) = -ig J dyMy)rG;(y)My) , (1) 

where ipg^gi are the wave functions of the initial quark and final quark, G is the wave 
function of the emitted gluon. We write each quark wave function in the form ijji{y) — 
e::fj^\—iEi{t—z)\u\(\)i{z., p) I \/1Ei (hereafter the bold vectors denote the transverse vectors), 
where A is quark helicity, u\ is the Dirac spinor operator. The z-dependence of the 
transverse wave functions 0j is governed by the two-dimensional Schrodinger equation 

« — Q-^ — = 1 -+9Qn[G^-GJ\(f>i{z,p), (2) 

where now G denotes the external vector potential (the superscripts are the Lorentz in- 
dexes and n = 1, 2 correspond to the A and B color components in the Cartan subalgebra), 
Qn is the quark color charge. The wave function of the emitted gluon can be represented 
in a similar way. We will assume that in the QGP for the usual parton masses one can 
use the corresponding quasiparticle masses. 

We take the external vector potential in the form Gf^ = [H„ x p]^, G„ = 0, G^ = 
(we assume that chromoelectric field is absent, however, it can be included as well). For 
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this choice of the vector potential the term —gQnG\ in (2) can be viewed as the potential 
energy in the impact parameter plane Ui = — Fj ■ p, where Fj is the corresponding Lorentz 
force. Then the solution of (2) can be taken in the form 



(pi{z,p) = exp\ipi{z)p 



2Ei 



Jo 



+ m„ 



(3) 



Here the transverse momentum Pi{z) is the solution to the parton equation of motion in 
the impact parameter plane 

"^^^ Uz). (4) 



dz 



Below we denote the value of Pj(±oo) as pf. By using (1), (3) one can obtain 



{gq'\S\q)^-ig{2nfS{Eg + E,, 
X exp < —I I dz' 



oo 



^pl{z')+m'i 



dzV{z, {X})S{pJz) + pJz) - pJz)) 



pl{z')+ml pl,{z')+ml 



2E, 2Eg 



2Eg/ 



(5) 



where V is the spin vertex factor, {A} is the set of the parton helicities. For the transition 
conserving quark helicity V = —iEq\/l — x[2XqX+{2 — x)\g] [vx{z) —i\gVy{z)]/ \/2, and for 



the spin-flip case V — ixmq{2\q\g+l)/ sj2{l — x). Here \{z) = Vg{z)—Vg{z) is the relative 
transverse velocity in the final gq' parton system which can be written as v{z) — q{z) / ji 
with q(^) = Pg{z){l — x) — Pqi{z)x, fi = Eqx{l — x) (here x is the longitudinal gluon 
fractional momentum). Due to the color charge conservation Fg = Fg + Fq/. For this 
reason, the argument of the second 5-function does not depend on z, and can be replaced 
by pj^ + pi — pj'. From (5) we obtain the gluon emission spectrum 



X exp 




dP 

dx 

pl{z)+ml 
2E„ 



(2vr) 



dp'l / dzidz2g{zi, Z2) 



pl{z)+m\ 



2En 



pl,{z) + ml 

2En' 



where the vertex factor reads (we recover the vertex color factor C) 

Cn 

= .F2w/ 5: ^*(^2, {A})y(^i, {A}) = ^iv(;.2)v(^i) + g, 

QI^qX[l - X) 



(6) 



(7) 



with gi = Cq!s(1 — x + x'^/2)/x and g2 = CagTrigX^ /2iJ,'^ (the two terms in (7) correspond 
to the non-flip and spin-flip processes). The color factor reads C = |Ay^jX*/2p, where i, / 
are the color indexes of the initial and flnal quarks, Xa is the color wave function of the 
emitted gluon. 



For a uniform external fleld we can write v(z2)v(zi) = [q — f t /4]/ fi , where q = q{z), 
- {zi + Z2)/2), T = Z2 — zi, and f = dq/dz = Fg{l — x) — Fq/x. The argument of the 



exponential function in (6) can be rewritten as 

$(t, q) 



+ q2)r 



2/i 



+ 



fV 

24/x 



(8) 



with = rriqx'^ + rn'^{l — x). After replacing in (6) the integration over p+ by the 
integration over q we obtain for the radiation rate per unit length 



dP 



dLdx (27r)^ 



dq dr 




1 J —oo 





^2 



+ 92 



exp [-i$(T, q)] . 



With the help of integration by parts one can rewrite (9) as 



(9) 



dP 



dq dr 




1 J —oo 





dLdx (27r)2 
After integrating over q (10) takes the form 
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exp [-i$(T, q)] 



(10) 



dP 
dLdx 



iji r°° dr 

27r 7-00 T 




-92 



exp < —I 



2fx ^ 24// 



(11) 



Note that in (9)- (11) it is assumed that r has a small negative imaginary part. One can 
easily show that in (11) the integral around the lower semicircle near the pole at r = 
plays the role of the f = subtraction term. Expressing the integrals along the real axis 
in terms of the Airy function Ai{z) — ^y^i^i/a (22:^/^/3) (here Ki/^ is the Bessel function) 
(11) can be written as 

dP n , r°° 

-AiiK) + b dyAi{y), (12) 

K J K 



where a — —2€ gi/fi, b = fig2 



dLdx 

e^9i/li, K 



V(/.¥)V3. 



Prom (11) one can obtain for the coherence length of the gluon emission Lc ~ 
min(Li,L2), where Li = 2/x/e^ and L2 = (24///f^)^/"^. From this estimate one can easily 
show that for radiation of a charged gluon the condition L^/ Rg^i ^ 1 (which is neces- 
sary for validity of our small angle approximation) is really satisfied in the quasiclassical 
regime when EgRg ^ ^ 1 and Eg ^ nig. In the interesting to us region Fg/ < (hereafter 
Fi — |Fj|) for emission of a neutral gluon the condition Lc/Rq,L <S 1 is also satisfied. 

Our spectrum for neutral gluons when Fg = for nig = agrees with the photon 
spectrum obtained in the quasiclassical operator approach [17]^, and similarly to the 
photon emission (12) gives dP/dLdx oc x^"^^'^ at a; <^ 1. The nonzero gluon mass leads 
to the Ter-Mikaehan suppression at x < {rrig/EqFq/y^'^. In this region the parameter k 
becomes larger than unity and the spectrum dP/dLdx oc exp [— 2mg/3£'gFg/x^]/x. For 
the charged gluons from (12) one can obtain in the massless limit dP/dLdx oc x"^^"^ at 
a; ^ 1. The Ter-Mikaelian mass effect suppresses the gluon spectrum at x <j nig/EqFg. 
In this region dP/dLdx oc exp [—2ni^/3xEqFg]/x^^'^ at x <^ 1. Note that in general if one 
neglects the parton masses (or in the strong field limit) the spectrum (12) takes a simple 
form 

dP Q;,Cr(2/3)[l-x + xV2](9f')V3 



dLdx 



7rV^x[Eqx{l - x)]V3 



(13) 



^Note that for the photon emissfon by an electron the integrand in the exponential term in (6) can 
also be written as E^k^ x ^^^{t) / {E^ — E^) , where fc'' is the photon four momentum and x^{t) is the classical 
electron trajectory. This gives precisely the spectrum in the form obtained by Baier and Katkov [17]. 



4 



Our derivation is valid for the g ^ gg transition as well. In this case a and h have the 
same form but with = m^gil — x + x^) , gi = Cas[l + x^ + (1 — a;)^]/4a;(l — and g2 — 0. 
The gluon vertex color factor reads C = \XaXbXcfabc\'^ , where the index a corresponds to 
the initial gluon, and b, c to the final gluons. 

3. In our formula for the spectrum for a ^ be transition all the properties of the final be 
state are only accumulated in f = Flxl — 2FbFcXbXc + F^xl (except for the trivial vertex 
factor). The crucially depends on the relation between the forces acting on the final 
partons. The value of characterizes the difference in bending of the trajectories of the 
b and c partons in the external field which is responsible for the synchrotron radiation. 
For this reason the spectrum vanishes if at some x i = 0. This, for example, occurs for 
the gx — gygz gluon process at a; = 0.5 for the external field in the color state A. 

Note, that the spectrum obtained in the present paper can also be derived making use 
the light-cone path integral formalism [13] for gluon emission due to multiple scattering. In 
[13] the spectrum was expressed in terms of the Green's function for the Schrodinger equa- 
tion describing the relative motion in the qq'g system. In the absence of the external field 
the corresponding Hamiltonian has a purely imaginary potential U{p) = —inaqq'g{\p\)/2 
(here agqig is the cross section for the qq'g state, and n is the number density of the 
medium) . In the case of the synchrotron radiation (without multiple scattering) this po- 
tential should be replaced by the real potential U{p) = — f ■ p. The advantage of the 
path integral approach is that it allows one to treat the gluon emission including both 
the multiple scattering and bending of the trajectories in the external field. This analysis 
will be given in further publications. 

Our formula (11) disagrees with that obtained by Shuryak and Zahed in the soft 
gluon limit within the Schwinger's proper time method. [18]. In the spectrum derived 
in [18] (Eq. (20) of [18]) the argument of the exponential contains (we use our notation) 

2 *? 2 2 2 

FgiXg + Fg lustcad of our f . Also, in the pre-exponential factor instead of f there appears 
Fg,x'^g. Thus, our result at -C 1 agrees with that of [18] only for the QED hke processes 
(emission of the neutral gluons). For a real QCD process, with charged gluon and quark, 
due to the absence of the interference term the spectrum of [18] is insensitive to the 
relation between the signs of the color charges of the final partons. It is strange enough, 
since the difference in the bending of the final parton trajectories (which is responsible 
for the synchrotron radiation) is sensitive to the relation between the color charges of the 
final partons. Also, Eq. (20) of [18] gives clearly wrong prediction that in the massless 
limit the spectrum of the qi — > gzQs transition for the chromomagnetic field in the color 
state A vanishes (since in this case F^/ = 0). Indeed, this process except for the spin 
effects is analogous to the synchrotron radiation in QED, and there is no physical reason 
why it should vanish. Note also that since in [18] the pre-exponential factor contains the 
Lorentz forces acting on the final partons in non-symmetric form it is clear that for the 
g ^ gg process the method of [18] should give the spectrum with incorrect permutation 
properties Thus, one sees that the formula obtained in [18] clearly leads to physically 
absurd predictions. Unfortunately, the details of the calculations have not been given in 
[18]. For this reason it is difficult to understand what is really wrong in the analysis [18]. 

4. We use the quasiparticle masses obtained from the analysis of the lattice data within 
the quasiparticle model [19]. For the plasma temperature T ~ (1 — 3)Tc which is relevant 
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to the RHIC and LHC conditions the analysis [19] g 0.3 and rUg ^ 0.4 GeV. 

In Figs. 1, 2 we present the averaged over the color states gluon spectra for q gq' and 
9 ^ 99 processes for the chromomagnetic field in the color state A for different initial 
parton energies. The computations are performed for = 0.3 and gHA/m% = 0.05, 0.25 
and 1, where is the Debye mass (we assume that as for an isotropic weakly coupled 
plasma mjj — ^m^). The results for the chromomagnetic field in the color state B are 
very close to that shown in Figs. 1, 2, and we do not show them. The decrease of the 
spectra at a: — (and x ^ 1 for 9 ^ 99 process) which is well seen for the smallest 
value of the field is due to the Ter-Mikaelian mass effect. This suppression decreases with 
increase of the chromomagnetic field. 

In Fig. 3 we show the spectra for the q — > gq' process for different quark and gluon 
color states. To illustrate better namely the dependence on the color indexes in all the 
spectra we replace the color vertex factor C by unity. One can see that the spectra grow 
strongly in the region of small x with the gluon color charge. For this reason the averaged 
over color spectra shown in Figs. 1, 2 are dominated by the processes with charged gluons. 

For understanding the potential role of the synchrotron- like radiation in the jet quench- 
ing in AA-coUisions it is interesting to estimate the energy loss due to the synchrotron 
gluon emission. Of course, a realistic estimate of this effect requires detailed information 
on the time evolution of the QGP instabilities. Also, as we said in the introduction, an 
accurate analysis should treat the synchrotron radiation and usual bremsstrahlung due 
to multiple scattering on an even footing. This is, however, beyond the scope of the 
present analysis. In the present paper we can give only a crude estimate of the effect. 
To estimate the chromomagnetic field we rely on the idea that the small viscosity of the 
QGP observed at RHIC is due to parton rescatterings in the turbulent magnetic field [10]. 
In this model the viscosity/entropy ratio r]/s ^ 1/g'^^^^'^ [10], where ^ is the anisotropy 
parameter of the initial plasma distribution. It is expected that the generated magnetic 
field is saturated at g'^{H'^) ~ C'^iTi'j^. To obtain a realistic r]/s ratio one should assume 
that ~ 1. Probably, the value of the magnetic field obtained in this way may be viewed 
as an upper bound. Indeed, in this case the ratio of the magnetic energy to the thermal 
parton energy is ~ 0.3, and a higher fraction of the magnetic energy looks unrealistic. 
Making use the estimated magnetic field for RHIC conditions for — 0.3 we obtained 
AE/E ~ 0.1 - 0.2 for quarks and AE/E ~ 0.2 - 0.4 for gluons at ^ 10 - 20 GeV 
(for as = 0.5 the results are about two times bigger). These estimates have been ob- 
tained assuming that the parton path length in the magnetic field is ~ 2 — 4 fm. Also, 
we neglected any finite-size effects. These effects may be important if Lc and L are of 
the same order. For RHIC conditions the dominating contribution to the energy loss 
comes from the soft gluon emission where Lc ~ 1 — 2 fm. In this situation the finite-size 
effects may suppress the energy loss by a factor ~ 0.5. One more suppression mechanism 
may be connected with the finite coherence length of the turbulent magnetic field, Lm- 
However, if for the unstable magnetic field modes the wave vector k"^ < ^mjy [10], for 
RHIC conditions this suppression should not be very strong since we have L^/Lf. ^ 1. In 
this regime the turbulent effects should not suppress strongly the energy loss, and as a 
plausible estimate one can take the turbulent suppression factor ~ 0.5. Even with these 
suppression factors the synchrotron energy loss turns out to be comparable with that due 
to ordinary bremsstrahlung, and its relative contribution may be larger than that from 
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the coUisional energy loss [20]. Of course, the above estimates are very crude. Neverthe- 
less they demonstrate that the synchrotron radiation can be important in jet quenching 
and deserves further more accurate investigations. Of course, the synchrotron radiation 
and bremsstrahlung are not additive, since magnetic bending of the parton trajectories 
will suppress bremsstrahlung due to multiple scattering and vice versa multiple scattering 
will suppress the synchrotron radiation. This, in principle, can mask the effect of the syn- 
chrotron emission in the energy loss. However, the synchrotron radiation can reveal itself 
in the longitudinal broadening of the jet cone. Note that for the synchrotron emission this 
effect appears already at the level of the gluon emission itself contrary to the mechanism 
of Ref. [11] where it is a purely final-state interaction effect. 

5. In summary, we have developed a quasi classical theory of the synchrotron-like gluon 
radiation. In the QGP the gluon spectrum is dominated by the processes with emission 
of the charged gluons, the effect of the neutral gluons is relatively small. Our calculations 

show that the parton energy loss due to the synchrotron radiation may be important in 
the jet quenching if the QGP instabilities generate magnetic field {H"^) ~ m'jj/g'^. Our 
gluon spectrum disagrees with that obtained by Shuryak and Zahed [18]. We give simple 
physical arguments that the spectrum derived in [18] is incorrect. 
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Figures 




Figure 1: The spectrum for the q gq process in the chromomagnetic field in the color 
state A for as = 0.3, gHA/mj^ = 0.05 (a), 0.25 (b) and 1 (c), for the initial quark energies 
Eg = 20 GeV (black), Eg = 40 GeV (red). Eg = 80 GeV (blue). 
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Figure 2: The same as in Fig. 1 but for the g gg process. 
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X 



Figure 3: The spectrum for the q ^ gq process for specific color states. The curves 
correspond to qi QAqi (black), qi gxq2 (red), qi gzqs (blue), qs gyqi (green). 
In all the processes the vertex color factor C in (7) is replaced by unity. The computations 
are performed for = 0.3, and Eg = 20 GeV. The chromomagnetic field is in the color 
state A with gHA/m% = 0.05 (a), 0.25 (b) and 1 (c). 
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